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SUMMARY 


A decoupled control approach for a Large Gap Magnetic Suspension System (LGMSS) is 
presented. The control approach is developed for an LGMSS which provides five-degree-of-freedom 
control of a cylindrical suspended element that contains a core composed of permanent magnet 
material. The suspended element is levitated above five electromagnets mounted in a planar array. 
Numerical results are obtained by using the parameters of the Large Angle Magnetic Suspension 
Test Fixture (LAMSTF) which is a small scale laboratory model LGMSS. 

INTRODUCTION 

This paper presents a decoupled, single-input single-output (SISO) control approach for a Large 
Gap Magnetic Suspension System (LGMSS). The control approach is developed for an LGMSS 
which provides five-degree-of-freedom control of a cylindrical suspended element that contains a 
core composed of permanent magnet material. The suspended element is levitated above five 
electromagnets mounted in a planar array. The LGMSS is a conceptual design for a ground-based 
experiment which could be used to investigate the technology issues associated with: magnetic 
suspension at large gaps, accurate suspended element control at large gaps, and accurate position 
sensing at large gaps (ref. 1). This technology would be applicable to future efforts which range 
from magnetic suspension of wind tunnel models to advanced spacecraft experiment isolation and 
pointing systems. An analytical model of an LGMSS configuration using five electromagnets 
mounted in a planar array is developed in reference 2, This model is used to investigate two 
LQR control approaches for the LGMSS in reference 3. In reference 3, the simplifying assumption 
is made that the change in field and field gradients with respect to suspended element displacements 
is negligible. In reference 4 the analytical model developed in reference 2 is linearized and extended 
to include the change in field and field gradients with respect to suspended element displacements 
and the open-loop characteristics of the resulting system are investigated. The purpose of this paper 
is to develop a decoupled SISO control approach for an LGMSS using the extended model developed 
in reference 4. The control approach is proportional-derivative (PD) where the command torques 
and forces are functions of positions and derivatives of position. Numerical results are obtained 
for a candidate design which uses parameters for a Large Angle Magnetic Suspension Test Fixture 
(LAMSTF). The LAMSTF was designed and built in order to investigate the feasibility of the 
LGMSS concept and to provide a test fixture for developing and demonstrating control approaches. 
Reference 5 presents a description of the LAMSTF and some of the control approaches which have 
been investigated. LAMSTF suspended element parameters and field components generated by the 
electromagnets at the centroid of the suspended element are presented in the Appendix. 
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SYMBOLS 

system matrix (state-space representation) 
input matrix (state-space representation) 

modified input matrix (eq. (13)) 
magnetic flux density vector 
total force vector on suspended element 
command force vector 
gain matrix (eq. (59)) 

forward-loop transfer function matrix (eq. (33)) 
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acceleration due to gravity (g « 9.81 m/s 2 ) 

suspension height, suspended element centroid to top plane of coils 
coil current vector 

suspended element transverse moment of inertia 

coefficient matrix of field or field gradient components 

bias stiffness terms defined by equations (24), (25), and (26) 
respectively 

bias stiffness terms defined by equation (22) 
bias stiffness terms defined by equation (23) 
position and rate gain for x control loop 
position and rate gain for y control loop 
position and rate gain for z control loop 

i 

position and rate gain for 9 y control loop 

position and rate gain for 9 Z control loop 

magnetization vector 

suspended element mass 

total torque vector on suspended element 

command torque vector 

velocity vector 

permanent magnet core volume 

weighting matrix (eq. (3)) 

coordinates in orthogonal axis system 

position commands for x, y, z control loops respectively 

damping ratio for x, y, z control loops respectively 

damping ratio for 9 y and 9 Z control loops respectively 

Euler orientation, 3, 2, 1 rotation sequence 

position commands for 9 y and 9 Z control loops respectively 

rate of rotation 

natural frequency of x, y, z control loops respectively 
natural frequency of 9 y and 9 Z control loops respectively 

Subscripts 


2 


electromagnet axes 

partial derivative of i component in j direction 
partial derivative of ij partial derivative in k direction 
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maximum value 
equilibrium condition 

components along x—, y—, z — axes respectively 
coil number 

Matrix Notation 

matrix 

inverse of matrix 

transpose of matrix 
column vector 
transpose of column vector 
row vector 


Dots over symbols denote derivatives with respect to time; a bar over a symbol indicates that it 
is referenced to suspended element coordinates. 


EQUATIONS OF MOTION 

The equations of motion for the LAMSTF axe developed in this section using the linearized, 
extended, open-loop model of an LGMSS configuration consisting of five electromagnets mounted 
in a planar array which is presented in reference 4. Figure 1 is a schematic representation of the 
LAMSTF configuration which shows the coordinate systems and initial alignment. The suspended 
element coordinate system consists of a set of orthogonal x, y, z body-fixed axes which defines the 
motion of the suspended element with respect to inertial space. The suspended element coordinate 
system is initially aligned with an orthogonal x, y, z system fixed in inertial space. A set of 
orthogonal x/,, y j, z^ axes, also fixed in inertial space, define the location of the electromagnet 
array with respect to the x, y, z system. The x^ and y& axes are parallel to the x and y axes 
respectively and the zj, and z axes are aligned. The centers of the two axis systems are separated 
by the distance h. The open-loop equations of motion axe of the form (ref. 2). 

X = /(X,u) (1) 

where the state vector X is 

X T — [Oy Oy Oz Vj* Vy Vz x y zj (2) 

and the input vector u is 

u T = |li h h h h\ (3) 

In reference 4 these equations are linearized around the nominal operating point X 0 , I 0 by 
performing a Taylor series expansion. The linearized equations are (neglecting second-order terms 

and subtracting out X 0 ) 

<5X = ASX + BSn (4) 
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where 


and 


A = vMxWi 




A 

= OX/dX 

1 







(5) 
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= dX/dl 
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where 


Wi = 


1 /le 0 

0 l/Ic 


l/m c 


l/m c 


l/m c 


(9) 


The K coefficients in equation (8) represent the magnitude of field or field gradient components, 
at X Q , produced by a given current denoted / max - These components are presented in Tables A2 
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and A3 in the Appendix. For further discussion see reference 4. From reference 4 equation (52), 


the equilibrium value of B xz is 

Bxz ~ m c g/vMx (10) 

Using the relationship of equation (10), element (5, 3) of the matrix in equation (7) reduces to 

m c g/vMx -2Bxz = -B xz (11) 

Substituting this value in equation (7) and expanding equation (4) results in 

Ic&y — vMxi-BxOy — B xz x ~ By Z y — B zz z) + (vM x /I max )(-[K z \{l}) (12) 

IcClz = vM x (— Bx@Z ”1" BxyX + Byyy + ByzZ ) + ( vM x / Tmax) ( [AyJ {I }) (13) 


mcVx — vM x (—B X z0y + 2 B xy 0 z + £( Xx )x x + -®(xx)j /V + B( xx } z z) + vM x j / ma x( L-^xxJ {7}) (14) 


m c Vy — vM x (By Z 6y + (Byy — B XX )6 Z + B^ X y^ x x + B( X y}yy + B^ X y^ z z) 

4 - vM x /I max {[_K X y\{l}) ( 15 ) 

tricVig — vM x ((B xx — B zz )9y + By Z 0 z + B( xz j x x 4 - B( xz ^yy + B^ xz ^ z z) 

+ vM x /I m ax{\_Kxz\{l}) ( 10 ) 


The first terms on the right in the above equations are the torques and forces generated on the 
core due to perturbations in X , evaluated in the presence of the uncontrolled fields and gradients 
produced by the constant or “bias” currents required to provide equilibrium suspension. The second 
terms are the torques and forces generated on the core by controlling the coil currents about the 
suspension currents. The controlled torques and forces can be written as 



(17) 


where 


[B] = (vM x )/I max 


Kzl 1 

L*»J 

lAzyJ 

LUCrjJ 


( 18 ) 


From equation (17), the currents required to produce given command torques or forces become 


{I> = [6r{y} 


(19) 
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This decouples the forces and torques in terms of commands. As can be seen in equations (12)— 
(16) however, the motion is still highly coupled through the bias terms. In order to determine the 
extent of this coupling, the values for the bias, or uncontrolled field and gradient terms, can be 
calculated using the suspension currents as discussed in reference 4. The suspension currents, from 
equation (54) in reference 4, are 


{lo} — Imax 


lKy\ 1 

-1 

° 1 



° 

[AxxJ 

4 

° 

[K X y\ 


° 

. \ Kxz\ - 


l B xz ) 


( 20 ) 


Using equation (10), with the suspended element parameters from the Appendix, the value of B xz 
required for suspension is found to be 0.0965 T/m. With this value of B xz , the currents required 
to provide equilibrium suspension were found to be (using eq. (20)) 


{/ 0 } T = [-14.169 -4.381 11.466 11.466 -4.381 J (Amp.) (21) 


These currents were then used to calculate the values of the bias fields and gradients. The results 
are presented in Table 1. Note that many of the terms are zero. Substituting these values into 
equations (12)— (16) and assuming that the control currents are generated by using equation (19) 


results in 

tly = vM x /I c {.m29 y - .0965a:) + (1 /I c )T yc (22) 

flf = vM x / J c (. 008202 ) + (1/ I C )T ZC (23) 

V x = vM x / m c (— .O9650 y + -4912a:) + (1 /m c )F xc (24) 

Vy = vM x /m c (M02y) + (l/m c )F y c (25) 

Vz = vM x /m c (—.009z) + (1 /m c )F zc (26) 


Since B( xx y z and B^ xz ^ x are very small, they were neglected in equations (24) and (26). For 
simplification define 


KB y o = KB Z0 = vM x {. 0082) = .0184 Nm (27) 

KB yex = KB xye = vM x (. 0965) = .217 N (28) 

KB X = uAfx(.4912) = 1.1048 N/m (29) 

KB y = vM x (MQ2) = 2.1146 N/m (30) 

KB Z = vM x {. 009) = .0202 N/m (31) 


Note that the positive bias terms in equations (22)-(26) cause open-loop instability. These terms 
are similar to the unstable bias flux stiffness terms encountered with small gap magnetic bearings 
which use permanent magnet bias flux (refs. 6 and 7). The positive bias terms in equations (22), 
(23), and (24) produce the unstable compass needle modes (modes 1 and 3) discussed in reference 4. 
Equations (22) and (24) are coupled through negative bias terms and this coupling produces the 
combined motion in modes 1 and 2. Mode 1, as mentioned, is unstable but mode 2 is a stable 
oscillatory mode. The positive bias term in equation (25) produces an unstable translation mode 
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(mode 5 in ref. 4). The negative bias term in equation (26) produces a stable oscillatory mode 
(mode 4 in ref. 4). The suspended element motions are uncoupled with the exception of pitch 
rotation ( 0 y ) and x translation (eqs. (22) and (24)). The open-loop eigenvalues are presented in 
Table 2. 


CONTROL SYSTEM EQUATIONS 

As mentioned earlier, the control approach is PD where the command torques and forces are 
functions of positions and derivatives of position. For PD control, the command torques and forces 
can be written as 


1/C — (KPy$ + sKRyff)(9yc — Oy) 

(32) 

n zc ^ (KP z ff + sKR z ,$)(0 z c “ @z) 

(33) 

Fxc = (KP X -t~ sKR x )(x c x) 

(34) 

Ffc = ( KPy + sKRy)(y c - y) 

(35) 

F-zc = {KP Z 4- sKR z )(z c — z) 

(36) 


where KP y $ , KP z q, KP x , KP y , and KP Z are position gains for 6 y , 6 Z , x, y, and z respectively, 
KR z 0 , KRyff, KR X , KRy, and KR Z are rate gains for 9 y , 6 Z , x, y, and z respectively, 
0 yc , 6 ZC , x c , I/O and z c are position commands, and s is the Laplace operator. A block diagram of 
this control system is presented in figure 2. 

Control of pitch rotation (0 y ) and x translation will be examined first since these are the only 
suspended element motions which are coupled. The approach is to close control loops around each 
axis independently and to determine the effect of the cross-coupling on the design of the individual 
loops. In matrix form, (32) and (34) can be written as 


ft Half's 


(Oy 


~ Oy) \ 

-x) I 


(37) 


where [Gp] is a forward-loop transfer function matrix defined as 

\Gf\ = 


{KP y 6 + sKRyO 0 

0 {KP x + sKR x )\ 


(38) 


Taking the Laplace transform of equations (22) and (24) and putting them in matrix form results 
/ A \ _ r (1// c)KB v , -(\/Ic)KBy0 I '\ ( Oy 1 f (l/lc)Tyc \ f39 , 

\i 2 ij [-(l/m.JJfBjto (1 /m c )KB x J \ x J \ (1 /m c )F ic } ( J 

Substituting for T yc and F xc in (39) and collecting terms results in 

Ul c J + KR y ,, + (KPy e -KB ve )) KB ** 1 { } = [G F ] { £} (40) 

^ K \ m c s Hr iv ti x & ~r *x A. & x )) J ^ x * a £ c ) 


The characteristic equation becomes 

(IcS 2 + KRyfS + (KPy,, - KByS )) 


det 


K ByQ X 


KByQ X 

(i m c s 2 + KR x s + (KP x - KB X )) 


= 0 


(41) 
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Expanding the determinant results in 

(a 2 + KR y 0s/! c * (KP yd ~ KB y6 )/I c )(s 2 + KR x s/m c + (KP X - KB x )/m c ) - KB 2 y6x /I c m c = 0 (42) 


The natural frequency of the Q y control loop, independent of the cross coupling, becomes 

u y0 = ~I C 

with a damping ratio of 

P* = KSyt/ii^KPy, - KBy,)I C ) 

Similarly for the x control loop 

w* = y/iKPx - KBx)JW c 
and 

p, = KR x /(2^/(KP x - KB x )m c ) 

The cross coupling term, KB^ X / I c m c , can be defined as an equivalent natural frequency 


(43) 

(44) 

(45) 

(46) 


ui c = KBygx/ii/lcmc) (47) 

Substituting (43)-(47) into (42) and expanding results in 

s 4 + (2 Py 0 U) y ff 4- 2p x u x )s z + (( 2 p y 0 Uyff)( 2 p x uf x ) + vlo+ u l) s 2 ( , 

+ ((2 Py0Uy0)(u%) + (2p x U x )(u}yg))s + U 2 0 U 2 - J 2 = 0 


The characteristic equations for the remaining, uncoupled, control loops can be obtained in a 
similar manner. For yaw rotation (9 Z ), y, and z translation respectively, the characteristic equations 
become 


s 2 + KR z0 s/I c + (. KP z e - KB z e)/Ic = 0 

(49) 

i s 2 + KRys/m c + ( KP y - KB y )/m c = 0 

(50) 

a 2 + KR z s/m c + (KP Z - KB z )fm c = 0 
Damping ratios and natural frequencies become 

(51) 

Pz8 - K R-z0/ ( 2 \j (KP zd - KB z e)l c ) 

(52) 

U, zS = yj(KPa - kbT 6 )/Ic 

(53) 

Py — KR y /{2 s j (KPy — KB y )m c ) 

(54) 

U)y = \J {KPy KBy)/m C 

(55) 

Pz = KR z /(2y/(KP z - KB z )m c ) 

(56) 

co z = yj (KP Z — KB z )/m c 

(57) 
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EXAMPLE DESIGN 

In order to investigate the effects of the cross-coupling discussed earlier, consider an example 
design where the natural frequency of each loop is set to 100 rad/sec and the damping ratio is set 
to 0.7. Using equations (43)-(47) and (52)-(57), the rate gains and position gains become 


KR yB — KR z0 = 7.7126(10) -4 Nms/rad 

(58) 

KR X = KRy = KR Z *= 3.0974 Ns/m 

(59) 

KPyff — KPzff = 0.0735 Nm/rad 

(60) 

KP X = 222.3447 N/m 

(61) 

KP y m 223.3545 N/m 

(62) 

KP Z = 221.2601 N/m 

(63) 


For simplicity, assume a zero set-point and perfect differentiators. The command torques and 
forces can then be written as 

r Ty C 'j KRy e 0 KP y0 0 0 0 0 0 0 0 ' 

Tic 0 KR zd 0 KR Z0 0 0 0 0 0 0 

1 F ic >= 0 0 0 0 KR X 0 0 KP X 0 0 {X} (64) 

Fyc J O 0 0 0 0 KRy 0 0 KPy 0 

,F- ZC J L 0 0 0 0 0 0 KR Z 0 0 KP Z 

or, using compact notation 

{£} = - [G|W (65) 

The current becomes, from equation (19) 

{/} = [BJ-'IGKX} (66) 

Substituting for u in equation (3) results in 

X = (A — BB~ l G)X (67) 

The eigenvalues of (A — BB~ l G) are presented in Table 2. The first two sets are related to the 
coupled axes (pitch rotation and x translation). The damping ratios and natural frequencies of 
these eigenvalues differ from the design values by approximately three percent. The same values 
are obtained try finding the roots of equation (43). Position and rate gains for pitch rotation 
and x translation can be independently varied to adjust the coupled damping ratios and natural 
frequencies. Varying these gains have no effect on the remaining eigenvalues. 

CONCLUDING REMARKS 

A decoupled, SISO, control approach for an LGMSS has been presented. The control approach 
is for an LGMSS which provides five-degree-of-freedom control of a cylindrical suspended element 
that contains a core composed of permanent magnet material. The control approach decouples the 
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five degrees of freedom in terms of force and torque commands. However, the pitch rotation and 
x translation axes remain coupled through bias terms resulting from the bias currents required to 
produce equilibrium suspension. Equations for position and rate gains for each axis are developed 
in terms of natural frequencies and damping ratios. In order to investigate the effects of pitch 
rotation and x translation coupling, an example design was performed using the parameters of 
the Large Angle Magnetic Suspension Test Fixture (LAMSTF) which is a small scale laboratory 
model LG MSS. Results of the design indicated that the damping ratios and natural frequencies of 
the coupled axes differed from design values but not by a significant amount (approximately three 
percent for the design considered). 
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APPENDIX 

LARGE ANGLE MAGNETIC SUSPENSION SYSTEM 
(LAMSTF) PARAMETERS 


This appendix presents, in the form of tables, LAMSTF suspended-element, permanent magnet 
core, and electromagnet parameters and components of fields and gradients (including second- 
order gradients) generated by the LAMSTF electromagnets at the centroid of the suspended 
element. The LAMSTF contains a planar array of five room-temperature electromagnets, with iron 
cores, mounted in a circular configuration. The configuration is shown schematically in figure 1. 
For a more detailed description of the LAMSTF see reference 5. The fields and gradients were 
calculated using VF/GFUN (ref. 8), including the pre- and post-processor OPERA, with all iron 
cores modelled. LAMSTF parameters are presented in Table Al. Electromagnet fields and first- 
order gradients are presented in Table A2 and second-order gradients are presented in Table A3. 
It should be noted that the full set of components is not included in the tables since B^ = 
and B^k ■— B^j^j = B^^. 


Table Al. LAMSTF Parameters 


Core diameter, m . . 8.509(10) -3 

Core length, m . . ................ 5.08(10) ~ 2 

Suspended element mass (m c ), kg . 22.124(10) -3 

Suspended element inertia (J c ), kg — m 2 . . 5.508(10) -6 

Core volume ( v ), m 3 ...... » 2.889(10)“ 6 

Core magnetization (M x ), A/m ...................... 7.785(10) 5 

Suspended element suspension height (h), m 0.1 

Electromagnet outer radius, m ..... 0.0825 

Electromagnet inner radius, m 0.0475 

Electromagnet height, m . . . 0.105 

Iron core radius, m . 0.038 

Location radius*, m 0.1375 


* Distance from center of array to axis of given coil. 
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Table A2. Electromagnet Fields and First-Order Gradients 
Calculated at Nominal Suspension Location 



B x 

By 

Bz 

Bxx 

B X y 

Bxz 

Byy 

By, 

B Z z 

Electromagnet 

(Tesla) 

(Tesla/meter) 

1 

+.0023 

0 

-.0009 

+.0218 

0 

-.0272 

-.0189 

0 

-.0029 

2 

+.0007 

+.0022 

-.0009 

-.0150 

+.0120 

-.0084 

+.0179 

-.0259 

-.0029 

3 

-.0019 

+.0014 

-.0009 

+.0077 

-.0194 

+.0220 

-.0049 

-.0160 

-.0029 

4 

-.0019 

-.0014 

-.0009 

+.0077 

+.0194 

+.0220 

-.0049 

+.0160 

-.0029 

5 

+.0007 

-.0022 

-.0009 

-.0150 

-.0120 

-.0084 

+.0179 

+.0259 

-.0029 


All values calculated for I max = 10 Amp. 


Table A3. Electromagnet Second-Order Gradients 
Calculated at Nominal Suspension location 



B(xx)x 

B(xy)x 

B(xz)x 

B( X y)y 

B[xy)z 

B(xz)z 

Electromagnet 



(Tesla/met 

er /meter) 



1 

+.0343 

0 

-.5347 

-.2192 

0 

+.0025 

2 

! 

-.1828 

-.1456 

+.1637 

-.1194 

+.0874 

+.0008 

3 

+.1656 

+.1373 

-.2679 

+.2290 

-.1413 

-.0021 

4 

+.1656 

-.1373 

-.2679 

+.2290 

+.1413 

-.0021 

5 

-.1828 

+.1456 

+.1637 

-.1194 

-.0874 

+.0008 


All values calculated for I max = 10 Amp. 
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Table 1. Bias Fields and Gradients 


Bxf Tesla . . . . , , . . , . 
B X x> Tesla/meter 

-0.0082 

0 

B X y> Tesla/meter ......... 

0 

I?ar*, Tesla/meter . . » , . . . . . 
Byy, Tesla/meter ......... 

By*, Tesla/meter . ........ 

B**, Tesla/meter 

0.0965 

0 

0 

. . . . . . . . ... . . . . . . . ... . . 0 

B(xx)xr Tesla/meter/meter .... 

B( XI )y, Tesla/meter/meter .... 

B( ia .) z , Tesla/meter/meter . , . . 

B( a .y) x , Tesla/meter/meter . . . . 

B(*y)y, Tesla/meter/meter . . . . 

B(a;y)*, Tesla/meter /meter .... 

B(xz)xi Tesla/meter/meter .... 

Tesla /meter /meter * . . . 

0.4912 

0 

0.0002 

0 

0.9402 

............ . . . 0 

0.0002 

0 

'*■1^ I fW* yf lAJ t * VUAI# 1 vv* 1 1AAV VV* ♦ • • * 

B(xz) z > Tesla/meter/meter .... 

-0.009 
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Table 2. System Eigenvalues for Example Design 


Open-loop 

Eigenvalues 

Damping Ratio ( p ) 

Natural Frequency (u n ) 

+/— 58.77 

Unstable 

58.77 

+/— 7.97i 

0 

7.97 

-h/— 57.8 

Unstable 

57.8 

+/— 0-96i 

0 

0.96 

+/— 9.78 

Unstable 

9.78 


Closed-loop 

Eigenvalues 

Damping Ratio (p) 

Natural Frequency (u/ n ) 

—70.01+/— 75.64i 

0.68 

103.1 

—70.01+/— 66.92i 

0.72 

96.8 

—70.01+/— 71. 42i 

0.7 

100 

—70.00+/— 71.43i 

0.7 

100 

—70.00+/ -71. 41i 

0.7 

100 
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